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The effect of the Hubbard interaction among conduction electrons on the double exchange 
model is investigated in a ferromagnetic metallic phase. Applying iterative perturbation theory 
to the Hubbard interaction within dynamical mean held theory, we calculate the one-particle 
spectral function and the optical conductivity, in which coherent-potential approximation is 
further used to treat the ferromagnetic Hund coupling between conduction electrons and localized 
spins. Identifying the decrease of the magnetization for the localized spin with the increase of 
the temperature, we discuss the temperature dependence of the one-particle spectrum and the 
optical conductivity. It is found that the interplay between the Hund coupling and the Hubbard 
interaction dramatically changes the spectral function, while it is somehow obscured in the 
optical conductivity. 
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§1. Introduction 

There has been a resurgence of interest in doped per- 
ovskite manganese oxides such as Ri_ K A x Mn03 (R is 
rare-earth and A is alkaline-earth ions) , since the discov- 
ery of the colossal magnetoresistance. These q 
show a variety of interesting phenomena 
which are caused by several competing interactions in- 



MM 



B b n t , al s o, char g e r aniL pi\ 
[00 B0 00 00001 



al de- 
Sincc 



eluding not onlv_ 
grees of freedomJlr 
the ferromagnetic state in a metallic phase is stabilized 
by the Hund coupling between localized spins and con- 
duction electrons, the double exchange model has been 
used as a proper model to describe the ferromagnetic 
metallic manganitesBl3# Various unusual properties ob- 
served experimentally have been clarified by using this 
model, in which the Coulomb interaction among; conduc- 
tion electrons has been neglected for simplicity ,□' lill EiP or 
has been assumed to be infinitely strongJi^ However, the 
finite Coulomb interaction should play an important role 
for the dynamics in Mn oxides, stimulating further the- 
oretical studies to take into account the finite Hubbard 
interaction. IL jshich have been done for the fully polar- 
ized case00B0EJ More recently, Held and Vollhardt 
have studied the effect of the finite Hubbard interaction 
on a paramagnetic metallic phase by means of dynamical 
mean field theory (DMFT) and have claimed that such 
correlation effects indeecLplay an important role together 
with the Hund coupling .Ef 

In this paper, we study electron correlation effects on 
the double exchange model in a ferromagnetic metal- 
lic phase by introducing the finite Hubbard interaction 
among conduction electrons. The.piiesent study is based 
on the approach of Shiba et al.ll3'L!J) with coherent po- 
tential approximation (CPA), and is in some respects 
complementary to the work done by Held and Vollhardt 
for a paramagnetic phase.ElP By exploiting iterative per- 



turbation theory (IPT) within DMFT for the Hubbard 
interaction, we compute the one-particle spectrum and 
the optical conductivity in which we further use CPA 
to deal with the Hund coupling. Based on the obtained 
results, we qualitatively discuss the temperature depen- 
dence of the above-mentioned dynamical quantities. It is 
shown that the introduction of the Hubbard interaction 
dramatically changes the one particle spectral function 
in the double exchange model. We also find that the 
optical conductivity is affected by the interplay of these 
two interactions, although it is not so conspicuous as in 
the one-particle spectrum. 

§2. Iterative Perturbation Theory Approach 

We study the double exchange model with the on-site 
Hubbard interaction in three dimensions. The Hamilto- 
nian is given by 



H 



J £' 



J7^n lT n a , (1) 



where c\ a {ci a ) is the creation (annihilation) operator of 
a conduction electron with the spin er(=T, [) at the site i, 
and iiia- — c\ a Ci a . For simplicity, the orbital degeneracy 
is neglected in this paper. The Hund coupling J between 
a localized spin Si (S = 1) and a conduction electron (<Tj 
is the Pauli matrix) is assumed to be ferromagnetic. 

We shall treat the electron correlations due to the Hub- 
bard interaction within DMFTpiJ which is known as a 
powerful method to systematically study strongly corre- 
lated electron systems. The treatment with DMFILbc.- 
comes exact in the limit of large spatial dimensions 
and even for three dimensions, it provides a powerful 
method which has successfully described a number of 
important phenomena, such as the Mptt transition, in 
strongly correlated electron systemsO'OEr To perform 
the DMFT procedure, we further make use of IPT (see 
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below),@0il'@0 which correctly reproduces the self- 
energy both in the weak coupling limit as well as the 
atomic limit, and furthermore properly interpolates the 
two limiting cases. 

We use CPA to treat the ferromagnetic Hund coupling. 
The validity of applying CPA to. the double exchange 
model is summarized as follows.EiP We are now dealing 
with the ferromagnetic Hund coupling, so that the dy- 
namics of localized spins may be described rather well 
by that of classical spins, because the ferromagnetic or- 
der considerably suppresses quantum effects. It may be 
thus legitimate to replace the local spin Si by the classi- 
cal Ising spin, and to further use CPA for treating Hund 
coupling for partially polarized ferromagnetic phase,tlP 
in which random potentials are induced by thermal fluc- 
tuations of the Ising spins. Note that the CPA procedure 
is properly embedded in the self-consistent equations for 
DMFTj£iP because both methods are based on local ap- 
proximations. 

We carry out the following self-consistent procedure 
for IPT combined with CPA. Since the original lattice 
problem is converted to an effective single impurity prob- 
lem in DMFT, let us first introduce the cavity Green 
function Q a (ui) and the dressed Green function G a {u>), 

g a (u) = [u + i6 + f i-E f(T -A a (Lo)]- 1 , (2) 

G <T H = [g- 1 H-S [/CT ( W )]- 1 , (3) 

where p and A a (u)) are the chemical potential and the 
renormalized hybridization function, respectively. Here 
we have denoted the self energy due to the Hubbard in- 
teraction as Ej/cr and Ef a — E f + Un s where Ef is the 
impurity level in the reduced single-site problem. Using 
the cavity Green function, the second-order self-energy 
due to the Hubbard interaction is expressed in a usual 
form, 



dei 
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„ /(ei)/(-e 2 )/(e 3 ) + /(-e!)/(e 2 )/(-e 3 ) 



(4) 



uj + id - ei + 62 — £3 

with Ptj(e) = — l/7rIm(C? CT (e)), where /(e) is the Fermi 
distribution function. For arbitrary filling, we improve 
the approximation by modifying the second-order self en- 



ergy as, 



with 



E[/ CT (w) 



(5) 



n a (l-n a ) (1 - n a )U + E f + p 

where riQ a and n a are the particle numbers respec- 
tively defined by = — I /it J di j/(ij)TrnG r (uj) and 
n a = -l/7r/^ o dw/(w)ImG(w)@0 ) Note that the co- 
efficients A„ and B a have been introduced so as to re- 
produce the correct results in the high energy limit and 
the atomic limit, respectively. 

We now incorporate the effect of the Hund coupling by 



introducing an additional self energy E# CT , which enters 
in the local Green function as, 



de- 



N (e) 



i5- e + p, - Ef a - H Ua (uj) - T, Ha (uj) 

(7) 

where G a (k,u>) denotes the one-particle Green function 
in the lattice system. We employ the semielliptic form 
as the bare density of states, No(e) = 2/ (irD 2 )^/ D 2 — e 2 
with the bandwidth D. Following a standard CPA proce- 
dure, the average of the random scattering of conduction 
electrons due to localized spins should be zero, leading 
to 

~JS l <7 - E gCT (o;) \ 

•1 - GJ? C M(- J Sid - E Hff (w)) Is.. 

which is rewritten in terms of the magnetization M of 
localized spins as, 



0, 



(8) 



GT{u)V Ha {u) + E ff(T (u/) 
- G 1 ° c (uj)(JS) 2 + oJM = 0. 



(9) 



The self-consistency condition requires the local Green 
function in Eq. (7) to be equivalent to the dressed im- 
purity Green function in Eq. (3). We thus obtain, 

a.m^^g^h + e^h. (10) 

Note that the chemical potential p should be determined 
by the Friedel sum rule in the reduced impurity problem. 
This completes the description of our IPT approach com- 
bined with CPA. 

§3. Correlation Effects on Dynamical Quantities 

We have performed the numerical calculation by iter- 
ating the procedure mentioned in the previous section 
until each quantity should converge within a desired ac- 
curacy. The one-particle spectrum of conduction elec- 
trons thus computed is shown in Fig. 1. We first mention 
that since we are dealing with the dynamics of localized 
spins by CPA, the change in the magnetization may be 
approximately considered to be caused by thermal fluc- 
tuations. El Namely, the magnetization of localized spins 
M is increased from (a) to (d), which implies in our treat- 
ment that the temperature is monotonically decreased 
from T = Tc (critical temperature) down to T = 0. 

We start our discussions with the case of (a), where 
there is no magnetization in the system. In this para- 
magnetic phase, the spectrum near the Fermi surface in 
the U = 2 case gets narrower than the non-interacting 
case. In addition to this band-narrowing effect, there 
is the life-time effect, which is accompanied by a small 
structure developed around the energy for the doubly oc- 
cupied state. These aspects are in. agreement with those 
obtained by Held and VollhardtEJ With the increase of 
M (decrease of the temperature), the spectral function 



Correlation Effects on the Double Exchange Model 



3 



shows quite different behaviors for up- and down-spin 
electrons. The spectral function for up-spin electrons 
is gradually reduced to the original semielliptic form, 
whereas that for down-spin electrons starts to develop 
several hump structures, whose origin will be discussed 
below momentarily. Note also in (b) that the spectrum 
with up- (down-) spin electrons in the U = 2 case shifts 
to a lower (higher) energy side compared with the nonin- 
teracting case, which is caused by the interplay between 
the Hubbard and the Hund interactions. Namely, since 
the magnetic susceptibility of conduction electrons is en- 
hanced by the Hubbard interaction, the magnetization of 
conduction electrons for a given M of localized spins be- 
comes larger for larger U via the Hund coupling, giving 
rise to the large self-energy shift observed in (b). We will 
see below that this shift affects the profile of the optical 
conductivity. For reference, we show the magnetization 
of conduction electrons as a function of M for several 
choices of U in Fig. 2. 

As M is further increased, the spectral function for 
down-spin electrons develops three distinct peaks (Fig. 
1(c)). To see the origin of the structure, we first recall 
that conduction electrons in (c) are almost fully polar- 
ized, though localized spins are still partially polarized 
(M = 0.8), as a consequence of the enhanced spin suscep- 
tibility (see Fig. 2). Thus adding a down-spin electron to 
an unoccupied site may have two possibilities depending 
on whether the corresponding localized spin is parallel or 
antiparallel to the spin of the added electron. The former 
contributes to the lowest-energy peak where the added 
electron gains the energy — J, while the latter forms the 
middle peak with the energy increment J for making an- 
tiparallel spins. On the other hand, if a down-spin elec- 
tron is added to the site already occupied by an up-spin 
electron, this makes the highest-energy peak in the spec- 
trum. Note that the position of the lowest-energy peak is 
pushed from the bare level up to slightly above the Fermi 
level, which is again caused by the self-energy shift due 
to the Hubbard interaction and the Hund-coupling. To 
clearly show how the above three-peak structure is de- 
veloped with the increase of the Hubbard interaction U, 
we display of M = 0.8 for several choices of U in Fig. 3. 
It is seen that the above-mentioned three-peak structure 
is gradually formed with the increase of U, which is ac- 
companied by the rather large self energy shift. Finally, 
in the case of (d), the localized spins are almost fully 
polarized, so that the spectrum for up-spin electrons is 
reduced to the non-interacting one, while the spectrum 
for down spin electrons still has the two-peak structure. 
As mentioned above, the lower peak represents a particle- 
addition spectrum to an empty site, while the upper peak 
corresponds to that for a singly occupied site. In all the 
cases (a)~(d), the interplay between the Hund coupling 
and the Hubbard interaction is important to determine 
the temperature dependence in the profile of the one- 
particle spectral function. 

Experimentally, several nontrivial effects on the pho- 
toemission spectrum haw been reported in the man- 
ganese compounds J2j'E£tE2P which may not be simply ex- 
plained without the Coulomb interaction. A typical ex- 
ample is that the band width observed in a paramagnetic 



phase is about twice as large as that expected from the 
band calculation. This may be caused by the correlation 
effects among conduction electrons, as recently explained 
by Held and Vollhardt with the use of DMFT.E2P Our 
results in Fig. 1(a) are consistent with theirs as well as 
the experimental tendency. Also, it has been reported 
that the coherent spectrum near the Fermi surface is in- 
creased as the temperature is decreased below the critical 
temperature, indicating the decreasejaLcffcctive electron 
correlations at low temperatures.oEj'ES* Such effects are 
indeed seen in the present results of Fig. 1 for electrons 
with up spin, namely, the decrease of the temperature 
gradually reduces the effective interaction among con- 
duction electrons in the filled band, being consistent with 
the experiments. Although our treatment here neglects 
several important effects due to the orbitally degenerate 
bands with anisotropic hoppings, etc, it turns out that 
some qualitative features obtained are consistent with 
those observed in the photoemission spectrum. 

In order to further investigate the dynamical response 
in the system, we calculate the optical conductivity. We 
employ the following formula for the optical conductivity 
in which the vertex correction is neglected, 

/OC />OC 
de' / de 
. - oo J — CO 

x N (e)A„(e, e')A a (e, lo + e') - /(f ' + ^ , (11) 
A a {e,uj) = -- 

IT 

x Im( r± Y 

^w + i(5-e + /i- Ef a - Sao-(w) - £#cr(cjy 

(12) 

Note that this expcessioiubecomes exact in the limit 
of large dimensions H3'I3eJ The computed results are 
shown in Fig. 4. As well known, the optical conductivity 
consists of two peaks, each of which corresponds to the 
Drude-like part and the interband part for which two sep- 
arated bands are formed by the Hund coupling. Compar- 
ing the U — 2 case with the non-interacting one, we can 
see that the structure near the u> ~ (Drude-like peak) 
gets somewhat narrower via the correlation effect due to 
the Hubbard interaction. Also, we find that the hump 
structure corresponding to the interband excitations are 
smeared by the life-time effect due to the Hubbard in- 
teraction U, and this effect becomes more conspicuous 
when U is large, being consistent with the results de- 
duced for a paramagnetic phase-EJ* It should be noticed 
here that the maximum position of the interband excita- 
tions in the U = 2 case shows a non-monotonic behavior 
with the increase of the magnetization M, which is not 
observed in the noninteractning case. Namely, the max- 
imum position in the U = 2 case is once lowered with 
the increase of the magnetization (M = 0.5), and then 
increased again up to the original position {M = 0.8). 
This is a consequence of the self energy shift induced 
by the Hubbard interaction and the Hund coupling: As 
mentioned in Fig. 1(b), the one-particle spectrum for 
U = 2 suffers from the large self-energy shift when M 
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is increased. If M is further increased, the self energy 
shift becomes small again, as seen in Fig. 1(c), (d), giv- 
ing rise to the non-monotonic behavior of the interband 
excitations as a function of M. If the Hubbard inter- 
action is further increased, the spectral weight of the 
interband excitations becomes small. In this way, the 
interplay of the Hubbard interaction and the Hund cou- 
pling appears in the optical conductivity, although it is 
somewhat obscured in contrast to a dramatic change in 
the one-particle spectrum. 

§4. Summary 

We have studied how the correlation effect among con- 
duction electrons affects dynamical properties of the dou- 
ble exchange model in a ferromagnetic metallic phase. 
In order to investigate dynamical quantities in the whole 
energy region, we have exploited iterative perturbation 
theory combined with coherent potential approximation 
for the Hund coupling. By approximately relating the 
decrease of the magnetization of the localized spin with 
the increase of the temperature, we have discussed the 
temperature dependence of the one-particle spectrum 
and the optical conductivity. It has been shown that 
the interplay between the Hund coupling and the Hub- 
bard interaction dramatically affects the profile of the 
one-particle spectrum, which results in a characteristic 
structure quite different from the non-interacting case. 
Some results have been found to be qualitatively con- 
sistent with photoemission experiments. On the other 
hand, the effect of the Hubbard interaction is somewhat 
smeared in the optical conductivity, although it indeed 
affects the profile for both of the Drude-like intraband 
part as well as the interband part. 

We have shown that there appears a non-monotonic 
temperature dependence in the one-particle spectrum 
as well as the optical conductivity in a ferromagnetic 
phase, which is caused by the self-energy shift due to 
the Hubbard interaction and the Hund coupling. For 
example, such a behavior emerges in the interband exci- 
tations in the optical conductivity in Fig. 4. As already 
pointed out, however, the interband peak-structure is 
smeared when the Hubbard interaction is increased, pre- 
sumably making it rather difficult to observe such a 
non-monotonic behavior. It seems thus interesting to 
check whether this type of the correlation effect could 
be observed experimentally in the one-particle density of 
states around the Fermi level. 

In this paper, we have neglected the orbital degeneracy 
for conduction electrons as well as the corresponding off- 
diagonal components for the hopping matrix. Correla- 
tion effects including the orbital degrees of freedom may 
be expected to provide further nontrivial properties. To 
this end, it is necessary to extend the present study to 
such orbitally degenerate cases in a ferromagnetic phase. 
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Fig. 1. Onc-particlc spectral function p(ui) calculated by IPT and 
CPA; (a) M = (Curie temperature), (b) M = 0.5, (c) M = 0.8, 
(d) M = 0.99 (T ~ OK). The hole (electron) concentration is 
x = 0.2 (0.8). In all figures, the left (right) panel represents the 
spectrum for up (down) spin electrons. The energy is measured 
from the Fermi level ui = 0, and the bandwidth D is taken to be 
unity. The Hund coupling is taken as J = 1, and the Hubbard 
interaction is set to be U = 2 (solid line) and U = (dashed 
line) . 
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Fig. 2. Plot of the magnetization for conduction electrons as a 
function of that for localized spins M for various U. The hole 
concentration is x = 0.2. The other parameters arc as in Fig. 1. 




Fig. 4. Optical conductivity a(ui) for (a) U = 0, (b) U = 2; M = 
(solid), M = 0.5 (dashed), M = 0.8 (dotted) and M = 0.99 
(dash-dotted), respectively. The insets show the low- frequency 
part. The other parameters are as in Fig. 1. 




